Dynamical trapping of light in modulated waveguide lattices 
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A discrete analogue of the dynamical (Kapitza) trapping eflfect, known for classical and quantum particles in 
rapidly oscillating potentials, is proposed for light waves in modulated graded-index waveguide lattices. As in 
the non-modulated waveguide lattice a graded-index potential can confine light at either normal or Bragg angle 
incidence, periodic modulation of the potential in the longitudinal direction enables to trap optical beams at 
both normal and Bragg incidence angles. © 2013 Optical Society of America 
OCIS codes: 190.6135, 230.7370, 020.1335 



Light propagation in photonic lattices has attracted a 
great interest over the past few years, with the observa- 
tion of a host of new phenomena with no counterpart in 
continuous media [1,2]. Flexible control of light trans- 
port and localization in such devices can be realized by 
breaking the translation invariance of the lattice along 
the propagation direction by either periodic axis bending 
or out-of-phase modulation of refractive index of adja- 
cent guides. Examples of light control in modulated lat- 
tices include dynamic localization [3-6], tunneling inhi- 
bition [7-9], multiband refraction control [10], polychro- 
matic diffraction management [11-13], and defect-free 
surface waves [14,15]. Most of such light control tech- 
niques bear interesting analogies with coherent control 
of driven quantum systems, such as electronic or matter 
wave transport in driven lattices [16]. 
In this Letter a mechanism of light trapping in graded- 
index modulated waveguide lattices is proposed, which is 
based on a discrete analogue of the Kapitza (or dynam- 
ical) stabilization effect of classical and quantum par- 
ticles in rapidly oscillating potentials [17-21]. Dynami- 
cal stabilization generally refers to the possibility for a 
particle to be trapped by a rapidly-oscillating potential 
in cases where the static potential cannot. Well-known 
paradigms are the Kapitza stabilization of the pendu- 
lum [17, 18] and Paul traps for charged particles [19]. 
Here it is shown that modulated lattices can trap light 
beams under conditions where the non-modulated lattice 
cannot. Let us consider light propagation in a modulated 
waveguide lattice [7], which in the tight-binding approx- 
imation is governed by the discrete Schrodinger equation 

^~tz ^ + V'n-l) + ^n{z)ll)n , (1) 

where tpn{z) is the light field amplitude trapped in the 
nth waveguide, z is the longitiidinal propagation coor- 
dinate, K > is the coupling constant between adjacent 
waveguides, and ^n{z) is the longitudinal modulation of 
the propagation constant for the nth guide. In the fol- 
lowing, we will assume <I'„(z) = f{z)Vn, where f{z) is 
a modulation function with spatial frequency cj = 2tt/A 
and zero mean, and Vn is the graded-index (static) po- 
tential. In waveguide arrays manufactured by femtosec- 



ond laser writing, the modulation can be realized by 
slightly varying the writing speed for each waveguide [7] . 
In the absence of the longitudinal modulation, i.e. for 
/ = 1, a graded-index potential Vn can be exploited to 
focus or trap light beams. Such a kind of focusing has 
been proposed, for example, to achieve deep subwave- 
length focusing in metal-dielectric waveguide arrays [22] . 
A schematic of the refractive index profile for a non- 
modulated graded-index array is depicted in Fig.l. To re- 
late the propagation properties of the graded-index mod- 
ulated lattice (1) with the dynamics of a quantum par- 
ticle in a rapidly-oscillating potential [20,21], it is worth 
observing that the solution to the coupled-mode equa- 
tions (1) can be written as tpn{z) = tp{x = n,z), where 
the continuous function tp{x,z) satisfies the Schrodinger 
equation i{d'ip/dz) = 'Htj} with Hamiltonian [23] 

H = -2kcos{p^) + f{z)V{x) (2) 

where Px = —id/dx and V{x = n) = Vn- For a po- 
tential V that varies slowly over one lattice period and 
for a wave packet with a narrow angular spectrum dis- 
tribution, the semiclassical (ray optics) equations for the 
mean values {x) of wave packet position (in units of the 
array period) and {px) of refraction angle, as given by 
the Ehrenfest theorem, read [23] 

f).. *J>.-/,.,(?^)(«,., 

(3) 

Let us assume that the static potential V{x) is a bell- 
shaped potential, for example, described by a parabolic 
or a Gaussian function, with a maximum at a; = 0, and 
let us first consider the non-modulated array (/ = 1). 
Then, according to Eqs.(3) there are two fixed points, 
((x) = 0, (px) = 0) and ((x) = 0, (p^;) = ti"), the for- 
mer being imstablc and the latter being stable [24] . This 
means that a broad beam launched into the array at 
nearly normal incidence will not be trapped by the the 
static potential Vm which acts as a defocusing lens for the 
discretized beam. Conversely, an injected broad beam 
tilted near the Bragg angle (corresponding to {px) — tt) 
will be confined by the potential y„, which acts as a fo- 
cusing lens. This is shown, as an example, in Figs. 1(a) 
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Fig. 1. (Color online) Beam propagation (snapshot of 
|'(/'„(z)p) in a modulated lattice (square-wave modula- 
tion, Lo / K = 0.5) with a Gaussian-shaped graded index 
potential (peak amplitude Vb) for normal and Bragg in- 
cidence angles and for increasing values of normalized 
potential peak amplitude: Vq/k = in (a) and (b) (non- 
modulated lattice); Vq/k = 3 in (c) and (d); Vo/k = 5 in 
(e) and (f). The right panels schematically show the in- 
dex profile of the non-modulated (upper plot) and mod- 
ulated (lower plot) graded-index waveguide arrays. 

and (b). If the sign of the potential Vn is reversed, 
the stability of the two fixed points is interchanged. 
Such a behavior is related to the well-known reversal of 
diffraction sign for discretized light at normal or Bragg 
beam incidence [1,25]. In fact, in the former case (nor- 
mal incidence) the Hamiltonian % is approximated as 
Ti. ~ — — 2K-|-y(a;), whereas in the latter case (Bragg 
angle incidence) one has Ti, ~ -\- 2k + V{x). Let us 
now consider the modulated lattice, and show that, sim- 
ilarly to the dynamical trapping of classical or quantum 
particles in rapidly-oscillating potentials [17,18,20,21], 
dynamical beam trapping in the lattice can be realized 
for both normal and Bragg beam incidence. Indeed, prop- 
agation of a broad beam, at either normal or Bragg inci- 
dence angles, is governed by the Schrodinger-type equa- 
tion with an oscillating potential, namely 



.dip 



f{z)V{x)yj 



(4) 



sign applies to nor- 
For a rapidly os- 
setting ^p{x, z) = 

and 



where the upper (lower) 
mal (Bragg) incidence, 
cillating potential, after 

</.(x, z) exp [-iV{x) (/; d^fiO - Jo d^iO^ ^ 
applying standard averaging methods [20, 21], at 
leading order the evolution of the slowly-varying 
envelope 4>{x, z) is described by a Schrodinger-type 
equation with an effective static potential Ve{x), namely 




(5) 

and the overbar denotes a spatial average over the 
oscillation cycle. For a bell-shaped potential V and 
considering normal beam incidence, the effective po- 
tential Ve{x) comprises two potential barriers, which 
can support metastable (resonance) states (see, for 
instance, [21]). Hence, though the static potential (i.e. 
for / = 1) can not trap light beams at normal incidence 
near a; = 0, in the modulated lattice this is possible 
owing to the dynamical (Kapitza) stabilization effect. 
For beam incidence at the Bragg angle, the sign of both 
diffraction and effective potential are reversed, and thus 
trapping is possible as well. It should be noted that the 
strength of the confining part of the effective potential 
is usually very small (it scales as ~ l/w^ [20, 21]), 
and the observation of such a dynamical stabilization 
effect for broad light beams in waveguide lattices might 
require extremely long propagation distances, which are 
not accessible with current waveguide array set-ups. 
Remarkably, we found that the dynamical trapping ef- 
fect persists even for relatively small spatial modulation 
frequencies w, of the order or smaller than the coupling 
constant k, making dynamical trapping observable 
with current waveguide array set-ups. As an example, 
Fig.l shows dynamical trapping of a broad Gaussian 
beam, at both normal and Bragg incidence angles, for a 
Gaussian-shaped potential Vn = Vq exp[— (n/w)^] and 
for a square-wave modulation function [f{z) = 1 in one 
semicycle, and f{z) = —1 in the other semicycle], as 
obtained by direct numerical simulations of Eqs.(l) for 
uj/k — 0.5, w = 18 and for increasing values of Vq/k- 
Initial condition is c„(0) = exp(— n^/25) in (a), (c) and 
(e) (normal incidence), and c„(0) = (— 1)" exp(— n^/25) 
in (b),(d) and (f) (Bragg incidence angle). For a typical 
coupling constant of k ~ 0.5 mm^^ [7], a propagation 
length of 40 in Fig.l corresponds to a physical length of 
~ 8 cm. Dynamical stabilization, observed in Figs.l(c-f), 
is related to the existence of metastable (resonance) 
states of the modulated lattice, which can be revealed by 
a direct computation of the quasienergy spectrum and 
Floquet eigenstates of the coupled-mode equations (1) 
with periodic coefficients. The quasi-energies (Floquet 
exponents) /i and corresponding Floquet eigenstates 
'ipil^\z) are defined as the solutions to Eqs.(l) of the form 
4^)(z) = w(f)(z)exp(-i^z) with ui''\z + A) = ui''\z) 
and —uj/2 < fj. < uj/2. In a truncated lattice, a 
metastable (resonance) state to Eqs.(l) corresponds 
to a Floquet state which is strongly localized near 
n — 0. For a modulation function f{z) satisfying the 
symmetry condition f{—z) — f{z), as the one used in 
the simulations of Fig.l, it can be readily shown that, if 
%l)^\z) — u"n\z) e:icp{—ifiz) is a Floquet eigenstate with 
quasienergy ^, then ipir'^\z) = (— l)"iil''''(— z) exp(i^z) 
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is also a Floquet eigenstate with quasicncrgy — /i. 
Therefore, resonance states appear in pairs. For the 
modulated lattice with parameter values used in the 
simulations of Fig.l, a direct computation of the Flo- 
quet exponents (assuming a truncated waveguide array 
comprising 2N + 1 = 161 waveguides, from n = —N to 
n = N) shows the existence of one pair of metastable 
states with quasienergies ±0. 0284k at Vq/k ~ 3, 
and of two pairs of metastable states (quasienergies 
~ ±0.0533k and ~ ±0. 1424k) at Vq/k = 5. As a general 
rule, for a fixed modulation frequency the number of 
resonance states increases as the amplitude Vq of the 
static potential is increased, which is in agreement with 
the prediction based on the cycle-averaged model: in 
that case the number of resonances sustained by the 
double-barrier effective static potential (5) increases 
as Vq is increased. The profiles of the two resonance 
states at Vq/k = 3, corresponding to the dynamical 
trapping of Figs. 1(c) and (d), are depicted in Fig. 2. 

In this case, the metastable state ipriT^'' with negative 
quasienergy is mainly excited by a normal incidence 
input beam, whereas the metastable state ipil^^ with 
positive quasienergy is the mainly excited state when a 
beam incident at the Bragg angle is considered. 
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Fig. 2. (Color online) Numerically-computed intensity 
distribution \ipl^''''\z)\'^ of the two metastable Floquet 
states, at the plane z — 0, for the modulated waveguide 
lattice of Fig.l with Vq/k — 3. The insets in the figure 
show the detailed behavior of the real (upper panels) 
and imaginary (lower panels) parts of the two Floquet 
states with negative (left plots) and positive (right plots) 
quasienergy. 

In conclusion, a discrete analogue of dynamical trap- 
ping of classical or quantum particles in rapidly oscillat- 
ing potentials has been proposed for light waves in pho- 
tonic lattices. Such an analogy could be exploited to trap 
light beams at different incidence angles, thus increas- 
ing the coupling efficiency from a broad angular spec- 
trum light source. However, as compared to ordinary fo- 
cusing and trapping in non-modulated graded-index lat- 
tices, dynamical trapping does not sustain truly guided 
modes, and requires a more complex refractive index 



management. Possible applications could be envisaged, 
for example, in light trapping at the sub-wavelength 
regime [22]. 
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